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Abstract 

In this paper, we consider the spectral stability of spatially periodic traveling wave 
solutions of the generalized Korteweg-de Vries equation to long-wavelength perturba- 
tions. Specifically, we extend the work of Bronski and Johnson by demonstrating that 
the homogenized system describing the mean behavior of a slow modulation (WKB) 
approximation of the solution correctly describes the linearized dispersion relation near 
zero frequency of the linearized equations about the background periodic wave. The 
latter has been shown by rigorous Evans function techniques to control the spectral 
stability near the origin, i.e. stability to slow modulations of the underlying solution. 
In particular, through our derivation of the WKB approximation we generalize the 
modulation expansion of Whitham for the KdV to a more general class of equations 
which admit periodic waves with nonzero mean. As a consequence, we will show that, 
assuming a particular non-degeneracy condition, spectral stability near the origin is 
equivalent with the local well-posedness of the Whitham system. 



1 Introduction 

In the area of nonlinear dispersive waves the question of stability is of fundamental im- 
portance as it determines those solutions which are most likely to be observed in physical 
applications. In this paper, we consider the stability of the spatially periodic traveling wave 
solutions of the generalized Korteweg-de Vries (gKdV) equation 



(1-1) u t = u xxx + f(u) x 
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where u is a scalar, x, t S R and / is a suitable nonlinearity. Such equations arise in a 
variety of applications. For example, the case f(u) = v? corresponds to the well known 
Korteweg-de Vries (KdV) equation which arises as a canonical model of weakly dispersive 
nonlinear wave propagation [T7J [22]. Moreover, the cases f{u) = (3u 3 , = ±1, corresponds 
to the modified KdV equation which arises as a model for large amplitude internal waves 
in a density stratified medium, as well as a model for Fernm-Pasta-Ulam lattices with 
bistable nonlinearity pQ [4]. In each of these two cases, the corresponding PDE can be 
realized as a compatibility condition for a particular Lax pair and hence the corresponding 
Cauchy problem can (in principle) be completely solved via the famous inverse scattering 
transforrrQ. However, there are a variety of applications in which equations of form (|l.ip 
arise which are not completely integrable and hence the inverse scattering transform can 
not be applied. For example, in plasma physics equations of the form (jl.ip arise with a 
wide variety of power-law nonlinearities depending on particular physical considerations 
[16] [18] [T9] . Thus, in order to accommodate as many applications as possible the methods 
employed in this paper will not rely on complete integrability of the PDE (jl.ip . Instead, 
we will make use of the integrability of the ordinary differential equation governing the 
traveling wave profiles: this ODE is always Hamiltonian, regardless of the integrability of 
the corresponding PDE. 

The stability of such solutions has received much attention as of recently: for example, 
see [5], [6], [7J, [8], [9], [10], [12], [13]. Here, we are interested in the spectral stability to 
long wavelength perturbations, i.e. to slow modulations of the underlying wave. There is 
a well developed (formal) physical theory for dealing with such problems which is known 
as Whitham modulation theory [27J |28j. This formal theory proceeds by rescaling the 
governing PDE via the change of variables (x, t) t— > (ex, et) then uses a WKB approximation 
of the solution and looks for a homogenized system which describes the mean behavior of 
the resulting approximation. Heuristically then, one may expect that a necessary condition 
for the stability of such solutions is the hyperbolicity- i.e., local well-posedness- of the 
resulting first order system of partial differential equations. In order to make this intuition 
rigorous, one must study in detail the spectrum of the linearized operator about a periodic 
solution in a neighborhood of the origin in the spectral plane, and compare the resulting 
low-frequency stability region to that predicted by hyperbolicity of the formal Whitham 
expansion. 

The first part of this analysis was recently conducted by Bronski and Johnson [6]. There, 
the authors studied the spectral stability of spatially periodic traveling wave solutions to 
(jl.ip to long-wavelength perturbations by using periodic Evans function techniques, i.e. 
Floquet theory. In particular, an index was derived whose sign, assuming a particular non- 
degeneracy condition holds, determines the local structure of the spectrum near the origin 
in the spectral plane. This index arises as the discriminant of a polynomial which encodes 
the leading order behavior of the Evans function near the origin for such perturbations: that 
is, it describes the linearized dispersion relation near zero frequency of the corresponding 

1 More precisely, the Cauchy problem for equation can be solved via the inverse scattering transform 
if and only if / is a cubic polynomial. 
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linearized equations about the background periodic wave. Notice that the origin in the 
spectral plane is distinguished by the fact that the tangent space to the manifold of traveling 
wave profiles can be explicitly computed, and the generalized null space of the corresponding 
linearized operator is spanned by a particular basis of this tangent space (again, assuming 
a particular non-degeneracy condition is met). In particular, Bronski and Johnson showed 
that if one assumes the conserved quantities of the PDE (jl.ip provide good local coordinates 
for near by periodic traveling waves, then the generalized periodic null space of the linearized 
operator is generated by taking specific variations in the parameters defining the periodic 
traveling waves. As a result, the linearized dispersion relation near zero frequency can 
be explicitly computed in terms of Jacobians of various maps from the traveling wave 
parameters to the conserved quantities of the PDE flow. These Jacobians were then shown 
to be explicitly computable in terms of moments of the underlying wave in the case of 
power-law nonlinearity f(u) = u p+1 , p £ N (see [7]). 

The purpose of this paper is to carry out the second part of the program, connecting the 
formal Whitham procedure with the rigorous results of [6]. To motivate our approach, recall 
that in [28] Whitham caries out the formal modulation approximation for the cnoidal wave 
solutions of the KdV equation, in which everything can be evaluated explicitly in terms of 
elliptic integrals. As a first step in his analysis, Whitham introduces a periodic potential 
function cj) with the requirement u = 4> x where u is a given solution of the KdV equation. 
In effect, the introduction of this potential allows for a Lagrangian formulation of the KdV 
which can can analyzed via the modulation theory presented in [28]. However, notice that 
by requiring the solution u to be of divergence form one is forcing that the solution have 
mean zero over one period. This of course presents no problem in the case of the KdV 
equation since all solutions can be made mean zero by Galilean invariance, but it does pose 
a serious problem when trying to extend Whitham's method to more general equations of 
form (jl.ip which admit periodic waves of non-zero mean; for example, the approximation 
would not be valid for the well known elliptic function solutions of either the focusing or 
defocusing modified KdV equation. 

In order to consider more general equations of the form (11. ip , then, it becomes imperative 
to find a way to side step this restriction to mean zero waves. To this end we follow the 
novel approach suggested by Serre [26], which is towork with the original variable u and 
augment the resulting WKB system with an additional conservation law associated with 
the Hamiltonian structure of the gKdV; see also the recent work of Zumbrun and Oh [20] 
in which this approach was used in the viscous conservation law setting. As we will see, 
this approach not only works in the KdV case originally considered by Whitham, but also 
for the general Hamiltonian gKdV. In particular, we will see that the linearized dispersion 
relation predicted by the linearized Whitham system correctly predicts (assuming particular 
non-degeneracy conditions are met) the modulational stability of the given periodic wave. 
To this end, we compare the linearized dispersion relation coming from the modulation 
approximation to that derived by Bronski and Johnson using Evans function techniques. 
Unlike the works of Serre and Zumbrun and Oh |26|. I20j. however, we show the equivalence 
by an essentially different method using direct computation and identities derived in [6] 
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as opposed to direct linear-algebraic comparisons. The remarkable fact that the resulting 
homogonized system can be explicitly computed in terms of the conserved quantities of the 
PDE flow and the parameters defining the periodic traveling waves seems to be very special 
to the case of the gKdV equation (jl.ip in which the ODE defining the traveling wave profiles 
is completely integrable. We expect that the approach proposed by Serre to show this 
equivalence should work in this case as well, by using the fact that additional conservation 
laws gives rise to additional first order variations, but the direct method employed here 
seems simpler in the present context. 

The results in this paper therefore extend the original work of Whitham concerning 
the modulational stability of cnoidal waves of the KdV to more general periodic waves 
of equations of form (jl.lj) . and also rigorously validate the predictions of the Whitham 
expansion as regards linearized stability about a periodic wave. However, as pointed out in 
|20[ I21j . this is only the begining of the story. The importance of the connection between 
Whitham expansion and linearized stability comes rather from the predictive power of the 
formal expansion in deriving more detailed information about behavior of solutions. See also 
the earlier work of Schneider [23} [241 [25] for reaction-diffusion and other equations. Thus, 
as in the related contexts [231 EH ESI [201 EI] , we expect the validation of this connection to 
lead to further important results about these and other systems. 

The outline of this paper is as follows. In section 2 we recall the basic properties of the 
periodic traveling wave solutions of (II. lh described in [6] and [13] . In particular, we discuss 
a parametrization of such solutions which will be utilized throughout the present work. 
In section 3, we outline the results of Bronski and Johnson concerning the low frequency 
asymptotics of the periodic Evans function and describe the connection to the modulational 
stability of the underlying periodic wave. In section 4, we carry out the slow modulation 
(WKB) approximation to derive the Whitham system for the gKdV equation. In particular, 
we connect this formal homogenization procedure to the low frequency analysis of Bronski 
and Johnson. We then close with a discussion and open problems. 

2 Periodic Solutions of the gKdV 

Throughout this paper, we are concerned with the periodic traveling wave solutions of the 
gKdV equation. To begin then, we recall the basic properties of the periodic traveling wave 
solutions of (jl.ip : for more details, see [6], [7], or [T3] . 

Such solutions are stationary solutions of (jl.ip in a moving coordinate frame of the form 
x + ct and whose profiles satisfy the traveling wave ordinary differential equation 

(2.1) u xxx + f(u) x - cu x = 0. 

This profile equation is Hamiltonian and hence can be reduced through two integrations to 
the nonlinear oscillator equation 

2 

(2.2) ^ = E + au+^u 2 - F(u) 
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where F' = /, F(0) = 0, and a and E are constants of integration. Thus, the existence of 
periodic orbits of (|2.ip can verified through simple phase plane analysis: a necessary and 
sufficient condition is for the effective potential energy V(u; a, c) = F(u) — au — |u 2 to have 
a local minimum. It follows that the periodic solutions of (|2.ip form, up to translatiorH, a 
three parameter family of traveling wave solutions of (II. ip which we can parameterize by 
the constants a, E, and c. In particular, on open sets in M 3 = (a, E, c) the solution to (|2.2p 
is periodic: the boundary of these open sets corresponds to the equilibrium solutions and 
solitary waves (homoclinic/heteroclinic orbits). 

In addition, we make the assumption that there exist simple roots u± of the equation 
E = V(u;a,c) which satisfy U- < u+, and that V(u;a,c) < E for u £ (u_,u + ). As a 
consequence, the roots u± are C 1 functions of the traveling wave parameters (a, E, c) and, 
without loss of generality, we can assume that u(0) = u_. It follows that the period of the 
corresponding periodic solution of (|2.ip can be expressed by the formula 

(2.3) T = T(a,E,c) = V2 ^ =\f /p^w— T 

where integration over T represents a complete integration from u_ to u + , and then back to 
U- again: notice however that the branch of the square root must be chosen appropriately in 
each direction. Alternatively, you could interpret T as a loop (Jordan curve) in the complex 
plane which encloses a bounded set containing both u- and n+. By a standard procedure, 
the above integral can be regularized at the square root branch points and hence represents 
a C 1 function of the traveling wave parameters: see [6] for details. 

Notice that in general, the gKdV equation admits three conserved quantities^] which can 
be represented as 



(2.4) M(a,E,c) = [ u(x)dx = ^ / 

Jo 2 Jr 



u du 



2 h ^/E-V{u;a,c) 



(2.5) P(a, E, c) = f T u(x) 2 dx = 

Jo 2 J T 



y/E-V(u;a,c) 



„ x f T ful „, A , V2 f E - V(u; a, c) - F(u) , 

H(a,E,c)= l-E- F (u))dx = ?-f K ' / K ' -du 

Jo \ 1 J 1 Jr v E - V (u; a, c) 



representing the mass, momentum, and Hamiltonian, respectively. As above, these integrals 
can be regularized at the square root branch points and hence represent C 1 functions of 
the traveling wave parameters. As seen in [6j and related papers, the gradients of the 
period, mass, and momentum plan a large role in the stability of the periodic traveling 
wave solutions of (11. ip . For notational simplicity then, we introduce the Poisson bracket 



2 The translation mode is simply inherited by the translation invariance of the governing PDE and can 
hence be modded out in our theory. 

3 In the case of the KdV or mKdV, the complete integrability of the PDE implies the existence of an 
infinite number of conservation laws. We will discuss this case more carefully in the following sections. 
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like notation 

for two-by-two Jacobians, and {/, g, h} Xt y ;Z for analogous three-by-three Jacobians. 

It should be pointed out that in [6] it was shown when E ^ gradients in the period 
can be interchanged for gradients of the conserved quantities via the relation 

(2.6) EV a: E,cT + aV a: E,cM + ~V a ,E,c + Va,E,cH = 

where V a ,E,c = {9 a ,dE,d c ). Thus, although all results of this paper will be expressed in 
terms of Jacobians involving the quantities T, M, and P, it is possible to re-express them 
completely in terms of conserved quantities of the gKdV flow. Such an interpretation seems 
possibly more natural and desired from a physical point of view. 

3 Evans Function Calculations 

Now, suppose u = u(-;a,E,c) represents a particular periodic traveling wave solution of 
(jl.ip . As we are interested in studying the spectral stability of this solution to localized 
perturbations, we must study the linearized equation about uq, namely 

d x C[u]v = -v t 

where v G L 2 (M) and the operator C[u] = —d x — f'(u) + c is a self adjoint periodic Hill 
operator. Taking the Laplace transform in time yields the linearized spectral problem 

(3.1) d x C[u]v = [xv 

where \x represents the Laplace frequency. We refer to the background solution u as being 
spectrally stable if the L 2 spectrum of the linear operator d x C[u] is confined to the imaginary 
axis: notice the spectrum is symmetric about the real and imaginary axis. 

The spectral analysis in L P (M), 1 < p < do, of differential operators with periodic 
coefficients is known as Floquet theory. In this theory, one rewrites the linearized spectral 
problem (|3.ip as a first order system of the form 

(3.2) Y x =U(x,li)Y 

and lets be a matrix solution satisfying the initial condition $(0, fi) = I for all 

/jfC, where I is the standard 3x3 identity matrix. The monodromy operator M(/i), or 
the period map, is then defined to be 

M(m) :=*(T,m). 

Notice that given any vector solution Y of the above first order system, the monodromy 
operator acts via the formula 



M.(fi)Y(x, p) = Y(x + T, p) 
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for all x G R and /j E C. In particular, it follows that the L 2 spectrum of the operator 
c^/^ii] is purely continuous and \i G spec(d x £[u]) if and only if 

det(M(/i) - e iK I) = 

for some k£1. Such an equality is equivalent with the existence of a nontrivial solution v 
of ([HD such that 

(3.3) u(a; + mT) = e imK v(x) 

for all x G R and m G Z, which is equivalent with (|3.1|) admitting a nontrivial uniformly 
bounded solution. 

Following Gardner (see [9] and [10J), we define the periodic Evans function for our 
problem to be 

(3.4) D(fx, k) = det (M(/x) - e iK l) , (//, k) G C x R. 

The complex constant e* K is called the Floquet multiplier and the constant k (which is not 
uniquely defined) is called the Floquet exponent In particular, by the above discussion 
we see fi G spec(d x C[u]) if and only if D(fj,,n) = for some n G R. Moreover, since the 
coefficient matrix H(/i, x) in (|3.2|) depends analytically on //, it follows that D is analytic 
in the variables \i and k and is real whenever \x is real. 

From a computational viewpoint, however, it is more convenient to define the periodic 
Evans function by choosing a basis /z) }?=i of the kerneHof the linear operator <9;,;£[u] — 

jj, and noticing from (|3.4[) that 

(3.5) D(h,k) = (det(cD,(0,/i)| j 3 =1 )) _1 det (^(T, - e^(0, /i)]^) 

This view is particularly convenient in our case: as a consequence of the integrability of 
the traveling wave ODE (|2.ip one can easily verify that the set of functions {u x , ue, u a } is 
linearly independent and satisfy 

C[u]u x = 0, C[u]ue = 0, and C[u]u a = —1. 

The first of these is a reflection of the translation invariance of (jl.ip . Thus, by taking 
variations of the underlying periodic solution u(-;a,E,c) of (|2.ip in the traveling wave 
parameters generates an explicit basis for the (formal) kernel of d x C[u]. In particular, 
one can easily verify that u x and a linear combination of u a and ue constitute bonafide 
T-periodic elements of the kernel of d x C[u], and hence D(fj,,0) = C(|/i| 2 ) for \fi\ -C 1. 
Continuing, it follows again from (|2.ip that d x Cu c = —u x and hence an appropriate linear 
combination of u c , u a , and ue lies in the first T-periodic Jordan chain in the translation 
direction and hence D(fi,0) = 0(|/i| 3 ) for \fi\ -C 1. Moreover, one can verify that if the 

4 By (|3.3p . the Floquet exponent encodes a class of admissible perturbations. 

5 Here, we are considering the formal operator d x C[u] with out any reference to boundary conditions. 
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Jacobian {T, M, P} a ,E,c is non-zero, then linear combinations of u x , u a , ue, and u c generate 
the entire generalized null space of 5 x £[n] and hence the Evans function should indeed be 
of order 0(|^| 3 ) near the origin. 

Using perturbation theory one can then analyze the way these solutions bifurcate from 
the (fi, k) = (0,0) state by using (|3.5f) with <&i(x, 0), &2(%i 0), and <&3(x,0) corresponding 
to u x , u a , and ue respectively. Moreover, it follows that the quantity 



9 ^ , ^ 



At=0 



corresponds to a linear combination of u c , u a , and ue, and that the first order variations in 
the u a and ue directions can be computed via variation of parameters. The details of this 
asymptotic calculation were carried out in Lemma 2 and Theorem 3 of [6], which we now 
recall as the main result of this section. 

Theorem 1. Assume that {T, M, P} a ,E,c is non-zero. Then in a neighborhood of (//, ft) = 
(0, 0) the periodic Evans function admits the following asymptotic expansion: 

D(jjl, k) = -^-{T, M, P} a ,E,c + "-^f ({T, P} E ,c + 2{M, P} a , E ) + m 3 + 0(| M | 4 + k 4 ). 

Notice that from Theorem [TJ spectral stability in a neighborhood of the origin is equiv- 
alent with the cubic polynomial 

(3.6) P{y) = -y 3 + | ({T, P} E , C + 2{M, P} a , E ) - ~{T, M, P} a , E , c 

having three real solutions, and hence modulational stability can be inferred from the dis- 
criminant of the polynomial P. Indeed, if this discriminant is positive then the polynomial 
P has three real roots and hence, using the Hamiltonian structure of the linearized operator 
d x C[u], the spectrum in a neighborhood of the origin consists of the a triple covering of the 
imaginary axis. If the discriminant is negative, however, then P must have a pair of complex 
roots yi 2 with non-zero imaginary part and hence the spectrum in a neighborhood of the 
origin consists of the imaginary axis with multiplicity one, along with two curves which are 
tangent to lines through the origin with angle arg^yi^) to the imaginary axis. In particular, 
the underlying periodic traveling wave is modulationally unstable in the latter case. This 
is the main result of [6] concerning the stability of periodic traveling wave solutions of (jl.ip 
to long wavelength perturbations. 

It should be pointed out that by standard abelian integral calculations, the Jacobians 
arising in Theorem [1] were shown in [7] to be explicitly computable in terms of the traveling 
wave parameters and moments of the underlying wave over a period in the case of a power 
law nonlinearity f(u) = u p+1 , p G N. For example, in the case of the KdV equation with 
f{u) = u 2 it was shown that the discriminant of P takes the form 

disc(Pf-)) = (a3 '° r3 + a ^ T2M + ai,2TM 2 + a , 3 M 3 ) 2 



disc(£ - V{-;a,c)) 
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where the constants ay 3_j represent nonlinear combinations of the traveling wave param- 
eters a, E, and c. In particular, it follows that the underlying periodic traveling wave is 
modulationally stable provided the equation 

E-V(y;a,c) = E + ay + ^y 2 - ^y 3 = 

has three real solutions, which is equivalent with the constants (a, E, c) corresponding to 
a periodic orbit of the traveling wave equation (|2,ip . Since every periodic traveling wave 
solution of the KdV can be represented in terms of the Jacobbi elliptic function cn(x;7), 
it follows that the cnoidal wave solutions of the KdV are always spectrally stable to long- 
wavelength perturbations as predicted by Whitham [28J. A similar representation holds 
in the case of the modified KdV equations with f(u) = =tu 3 and hence a given periodic 
solution u(-;a,E,c) of the modified KdV is modulationally stable provided the polynomial 
equation 

E + ay + ^y 2 - = 

has four real solutions. While the same procedure can be implemented for other power law 
nonlinearities f(u) = u p+1 with p£N, one must resort to numerical studies to compute the 
desired moments of the underlying periodic solution. 

Finally, it should also be noticed that taking k = in the low frequency expansion of 
the Evans function provided in Theorem [T] suggests that the stability index 

sgn (D(/i,0)D(A,0)) 

for real numbers 0< / u^;l^;A<oois determined^] precisely by the Jacobian {T, M, P} a ,E,< 
As a result, the underlying periodic traveling wave solution is spectrally unstable to co- 
periodic perturbations if {T, M, P} a ,E,c < 0, and (see [6]) is spectrally stable to such per- 
turbations provided that {T, M, P} a ,E,c and Te > 0. In particular, the low-frequency ex- 
pansion in Theorem Q] provides no information for the stability to co-periodic perturbations 
in the case where the Jacobian {T, M, P} a ,E,c is zero, i.e. when the map 

(a, E, c) G M 3 -> (T(a, E, c),M(a, E, c), P(a, E, c)) 

is a local diffeomorphism. Thus, by (12. 6ft . when E ^ the nonvanishing of {T, M, P} a ,E,c 
is equivalent with the statement that the conserved quantities of the PDE flow defined by 
(jl.ip provide good local coordinates for nearby periodic traveling waves. This nondegen- 
eracy condition has been seen in the stability analysis of periodic gKdV waves in several 
other contexts: for example, it appears in the nonlinear (orbital) stability analysis of such 
solutions [7] [13] to periodic perturbations, as well in the instability analysis to transverse 
perturbations in higher dimensional models [15J. In particular, using the elliptic integral 
methods [7] it has been seen that the quantity {T, M, P} a ,E,c is generically nonzero for a 
wide variety of nonlinearities, including the most physically relevant examples of the KdV 



6 It was shown in [6] by elementary asymptotic ODE theory that D(A, 0) < for A > sufficiently large. 
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and modified KdV equations. As we will see in the next section, at the level of the Whitham 
system this condition is necessary for the resulting first order system of conservations to be 
of evolutionary type. 

4 Slow Modulation Approximation 

We now complement the rigorous results of the previous section with a formal Whitham 
theory calculation. In particular, we want to show that the linearized dispersion relation 
P(-) in (13, 6h can be derived through a slow modulation (WKB) approximation, and hence 
that the formal homogenization procedures suggested by Whitham [57] [25J and Serre [2B] 
correctly describe the stability of the periodic traveling wave solutions of (jl.ip to long 
wavelength perturbations. To this end, recall from Section [2] that the gKdV admits a four 
parameter family T of periodic traveling wave solutions for some triple (a,E,c). We can 
thus form the quotient space V :=T/TZ under the relation 

(uKv) <=^ (3£ €R;t; = «(•-£)) • 

and we have the class functions 

T = T(u), a = A(u), E = E(u), c = C(u). 

Similarly, since the conserved quantities are translation invariant it follows that M and P 
can be interpreted as class functions on V themselves. 

Let {oq,Eq,cq) correspond to a particular non-constant periodic solution u$ of (|2.ip in 
V, that is 

Uo(x) = v(x;a ,E ,co) 

for some v £ V, where ao = A(uq), Eq = E(uq), and Co = C(uq). In particular, from the 
discussion in Section [2] it follows that since u$ is nonconstant the projection 

u i — ^ u 

is locally a fibration (where the fibers are circles), and hence V is locally of dimension three. 
Now, consider equation (jl.lj) in the moving coordinate frame x + cot: 

(4.1) ut = u xxx + f(u) x - c u x . 

A periodic solution of the corresponding traveling wave ODE is a stationary solution of 
(|4.ip . i.e. is a solution with wave speed s = in this moving coordinate frame. Letting 
e > 0, we now rescale (I4.ip by the change of variables (x,t) \— > (ex, et) so that the rescaled 
equation takes the form 



(4.2) 



ut = £ 2 u xxx + f(u) x - c u x , 
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where, with a slight abuse of notation, x and t now refer to the slow variables ex and et. 

Following Whitham then [27\ [28J, we consider a solution with a WKB type approxima- 
tion of the form 

(4.3) u £ (x, t) = u° fx, t, + eu l (x, t, + 0(e 2 ) 

where y t— > u°(x,t,y) is an unknown 1-periodic function. It follows that the local period 
of oscillation of the function u° is e/d x (j), where we assume the unknown phase a priori 
satisfies d x (j) / 0. We now plug (|4.3p into (|4.2p and collect like powers of e. The lowest 
power of e present is e~ 1 which leads to the equation 

fadyV? = {4> x d y ) 3 u° + (4> x d y ) f (u) - cq {4> x dy) u°. 

By defining 

(4.4) s = -^r, and u = <f> x , 

<Px 

it follows that the 0(e _1 ) equation is precisely the traveling wave ODE (|2.ip in the variable 
uy with wavespeed cq — s; In a similar way as above we may define class functions S and £1 
on V, respectively, associated to the quantities s and oj. Moreover, we may choose u°(x, t, •) 
to be spatially periodic and satisfy the nonlinear oscillator equation 

(4.5) = E + au° + (u ) 2 - F(u°) 

where a = a(x,t), E = E(x,t), and s = s(x,t) are independent of y. That is, we identify 
for each (x, t) £ M 2 the projection of the function y ^ u°(x, t, y) into V as being a periodic 
traveling wave solution of the gKdV in of form 

n°(x, t, y) = u{y; A(u°),E(u°), c - S(u )) 

for some function u G V, where we have used the notation from Section [2j In particu- 
lar, notice that the parameters (A(u°), E(u°), S(ii )) defining vP in V depend on the slow 
variables x and t introduced above. 

Continuing, we find that the O(e ) equation reads 

d t u° = d x f(u°) - cd x u° + d x {cP 2 x d 2 yU°) +d y (---) 

Averaging over a single period in y and rescaling we find 

(4.6) d t (M(u°)n(u )) - d x G(u°) = 
where, using T = T(u°) for convenience, 

G(u°) = ±£ {f(u°) - cu° + d 2 y n°) dy. 
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Using (|4.5[) and the periodicity of u° in the variable y, it follows that 



G(u°) = i I (A(u°) - S(u°)u°) dy = (A- SMtt) (u°) 



Jo 



Moreover, the Schwarz identity (j) x t = (fitx provides us with the additional conservation law 



Now, as noted above, the manifold V has dimension three and hence equations (|4.6p and 
(|4.7I) do not form a closed system for the unknown function vP. In [28] this problem was 
remedied in the case of the KdV equation by restricting to mean zero periodic waves, which 
forces the equation count to be correct. However, as noted in the introduction, this require- 
ment can not be implemented in the general case considered here since the gKdV in general 
admits periodic waves of non-zero mean. In order to close the system then, we follow the 
approach of Serre [26] and augment the above equations with an additional conservation law 
arising from the Hamiltonian structure of (|1.1|> . The choice of this extra conservation law 
seems somewhat arbitrary where there is more than one to choose; for example, in the case 
of the modified KdV where complete integrability implies the existence of infinitely many 
conservation laws. To identify a useful conserved quantity in our calculations, notice from 
the Section [3] that the conserved quantity P must play a significant role in the modulational 
stability of the underlying periodic solution. With this as motivation, we notice from (14. ip 
that 



Substituting (|4.3|) into (|4.8p as before and collecting like powers of e, we find that the 0(e^ 1 ) 
equation is the traveling wave ODE equation for u , as before, multiplied by the profile u , 
and averaging the O(e ) equation over a single period in y yields the conservation law 



(4.7) 








(4.9) 




where 



(4.10) 




+ u°«° - F(u°) 



^K) 2 -|(n°) 2 ) dy. 



Again, using (|4.5p one finds that 



v?f(u°) - F{u°) 



co - s 



2 



K) 2 + 




( 



- E — u°ul 
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and hence (|4.10p reduces to 

Q(*°) = ~£ (&. (u^ + E(u )^ dy 

(4.11) = -(SPn + 2E){u°). 

The homogenized system (|4.6p . (|4,7p . and (|4,9p is a closed system of three conservation 
laws in three unknowns, since tt° belongs to V. In particular, (|4.6p - (|4.9p describe the mean 
behavior of the WKB approximation (|4.3p . This system is precisely the Whitham system 
we were seeking, and it can be written in closed form as 

(4.12) d t (MO, Pn, fi) (it ) — d x (A — SMQ, -SPQ - 2E, -SQ) (u°) = 0. 

We now wish to make a connection between the system (|4.12p and the modulational stability 
of the original solution no- To this end, we make the assumption that the matrix 

(4.13) ^ o) ' 



d(a,E,a)K*E,s)=(ao,Eofl) 

is invertible, which implies that nearby periodic waves in V can be coordinatized by the 
traveling wave parameters (a, E, c) near (ao, Eq, cq). In particular, we can parameterize the 
Whitham system in terms of the variables (a, E, s) near (ao, Eq, 0) and hence 

(4.14) d t (Moj, Pu, u) (a, E, s) -d x (a- sMuj, -sPlo - 2E, -suj) (a, E, s) = 0. 

where now we have dismissed the use of the class functions on V previously defined and 
work directly with the functions M, P, and T = lo^ 1 defined on M 3 defined in section [2l 

We now linearize (|4.14p about the constant solution (a, E, s) = (ao, Eq, 0) corresponding 
to the original solution uq chosen above. To begin, notice from the integral representations 
CEl- ([23]) that 

d s (M,P,cj} = -d c (M,P,u), 

and hence assuming the quantity 

, fd(Moj,Puj,uj) l \ 1 rm „, , „ 

d6t V d(a,E,s) ka,E,s)=(ao,E ,o) ) = ^{T, M , P} a , E , c (a , E , c ) 

is nonzero, the corresponding linearized system is of evolutionary type, i.e. the system (|4.12p 
is of evolutionary type provided that the conserved quantities of the flow defined by (|l.ip 
provide good local coordinates for the nearby periodic traveling waves. Recall from Section 
[2]that this property is generically true for several physically relevant nonlinearities, and is a 
major technical assumption in much of the current work on the stability of such solutions: 
see again [6], [7], |13j . and [15]. Moreover the Cauchy problem for the linearized system is 
locally well posed provided that it is hyperbolic, which is equivalent with the characteristic 
equation 

,„ „„- , / d(Mto,Puj,uj) indfa- sMlo,-sPlo -2E,-sw)\ . n „. „ 

4.15 det [n \ ' ~7f— L w s 1 K40 =0 

\ o{a,E,s) T o(a,E,s) J 
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having three real roots in the variable j^, i.e. the matrix 

8(Muj,Plo,lo) 



d(a- sMu, -sPuj - 2E, -su) 



a ,E ,0) 



should be uniformly diagonalizable with real eigenvalues. 

To make the connection between the above (formal) modulation theory calculation and 
the rigorous Evans function theory of the previous section, notice that (|4.15p is equivalent 
to the vanishing of 

~ { d(Mu,Pu,u) iKd(a-aMu,-aPu-2E,-8u) 

(4.16) P{fi, k; a, E, 9) = deb [ fi- 



d(a,E,s) T d{a,E,s) 

at the point (ao, Eq,0). is nonsingular. Moreover, a straightforward computation shows 
that 

/ / M a uj M e uj -M c uj 

P([i,K;a ,E ,0) = /i 3 det P a to P E u -P c u 

\\ -T a u? -T E ^ T c u? 

2fj? ( n ( in 



P [ -- ) + (2M a -P E - 2T C 




where P(-) is the cubic polynomial given in (|3.6p encoding the modulational stability of the 
periodic traveling wave solution it(-; ao, Eq, Co). At first sight, this result is disturbing as it 
suggests that the formal WKB/homogenization calculation does not accurately describe the 
modulational stability of such solutions. However, this apparent discrepancy can be easily 
explained using the integral representations of the period, mass, and momentum given in 
(|2.3l) - (12.5p . Indeed, notice from this representation that 



M a = P E = 2T C 



V2 f v? du 



4 J r (E-V(u;a,c)) 3/2 
and hence we have the relation P E + 2T C — 2M a = 0. In particular, it follows that 

(4.17) P(v,K;a o ,E o ,0) = -^p(-j 

Notice that the root of the unexpected fj, 1— > — fx conversion in (I4.17h stems from the fact 
that in the Evans function calculation described in Section [3] we considered left moving 
waves (waves constant in the moving coordinate frame x + ct), while in homogenization 
calculations of the above form it is customary to use right moving waves. Indeed, from our 
definition of s in (14. 4p it follows that the quantity <j>{x — st) is constant, i.e. 4> corresponds 
to a right moving plane wave. Since changing right moving waves to left moving waves is 
equivalent to time reversal, which in turn arises as a/jn> — jj, at the level of the linearized 
equations, this added negative sign seems necessary in the framework presented here. In 
particular, we have proven the following theorem. 
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Theorem 2. Assume that {T, M, P} a ,E,c is non-zero and that the matrix in (|4.13p is non- 
singular. Then there exists a nonzero constant Tq such that the periodic Evans function 
admits the following asymptotic expansion in a neighborhood of ([i,k) = (0,0): 

D(fi, k) = T P(-fjt, k; oo, E , 0) + 0{\^ + k 4 ). 

That is, the linearized dispersion relation arising from the formal WKP approximation 
correctly describes the true linearized dispersion relation of the spectrum of the linearization 
about the underlying periodic traveling wave at the origin. 

Corollary 1. Under the assumptions of Theorem^ a necessary condition for the spectral 
stability of the periodic traveling wave solution uq(-) = u(-;ao,Eo,co) is that the spectrum 
of the three-by-three matrix 

„ s ( d{a- sMuj,-sPlo -2E,-suj)\~ l ( d{Muj,PLO,Lo)\ 

A ^ E ^--={ ) { d(a,E,s) ) 

d{Muj,Pu),uj) 

~~ d (a - sMu, -sPuo - 2E, -suj) 

be real at (ao, Eq,0): equivalently, that the Whitham system (|4.14p be hyperbolic at (ao, Eq,0). 

Proof. First, a straight forward calculation shows that the matrix 

d (a + sMuj, sPuj - 2E, -suj) 



d(a, E, s) 



(a,E,s)=(a ,E ,0) 



is invertible provided {T, M, P} a ,E,c is non-zero at (a,E,c) = (ao,E ,co) and the matrix 
in (|4.13j) is nonsingular at this special point. Moreover, if the matrix A(ao, Eq,0) had a 
nonzero eigenvalue r], then it follows from Theorem [2] that there is a branch of spectrum 
bifurcating from the origin admitting and asymptotic expansion 

» = ^ + 0(J) 

for | | <C 1. Thus, if r] has non-zero imaginary part one immediately has a modulational 
instability of the underlying periodic traveling wave solution. □ 

It follows that we have established that a necessary condition for the spectral stability 
of a periodic traveling wave solution of (|l.ip is that the Whitham system (|4.6p . (|4.7p . 
and (14. 9p be hyperbolic at the corresponding solution. Moreover, from the analysis in [6], 
if the matrix ^4(ao, -Eo> 0) has only simple real eigenvalues, it follows that the underlying 
periodic wave is modulationally stable. Indeed, the Hamiltonian structure of the linearized 
operator <9 x X[u] implies the spectrum is symmetric about the imaginary and real axes, and 
hence if the matrix A(a$, Eq,0) has real distinct eigenvalues a\, 02, and 03, then there 
are three branches fJ-j(n) which bifurcate from the origin with leading order expansions 
fij(n) = ^rf^ + 0(k 2 ). If one did not have spectral stability near the origin, then two of the 



5 DISCUSSION AND OPEN PROBLEMS 



16 



branches, say m and Hz, would satisfy Hi(k) = —^(k). It follows by equating like powers 
of k that a\ = az2, contradicting the simplicity of the spectrum of A(ao, Eo, 0). Thus, the 
strict hyperbolicity of the Whitham system is sufficient to conclude modulational stability 
of the underlying wave. 

5 Discussion and open problems 

In this paper, we have rigorously verified that the formal homogenization procedure intro- 
duced by Whitham to study the modulational stability of a periodic traveling wave profile 
of the gKdV equation does indeed describe the spectral stability near the origin. This in 
particular applies to the KdV case for which the WKB expansion is carried out in Chapter 
14 of Whitham [28] by Lagrangian methods making use of integrability to carry out com- 
putations in terms of elliptic integrals. (The rigorous connection between the Lagrangian 
formulation and the type of WKB expansion carried out here is an important result of |28j.) 

Recall from Section [J] that there is at first sight a "missing equation" in the WKB equa- 
tion for (KdV) and (gKdV). This is remedied in [28 J by introducing a "potential" consisting 
of the anti-derivative of the solution, thus viewing the equation as a one higher order PDE 
for which the equation count is correct. Here, we follow instead the approach suggested by 
Serre [26J of augmenting the system by an additional conservation law associated with the 
integral of motion coming from Hamiltonian structure. This approach works not only in the 
integrable KdV and mKdV cases, but also for the general Hamiltonian (gKdV). The deriva- 
tion of the Whitham equations seems also slightly simpler from this point of view. Indeed, 
using this approach of Serre together with identities developed by Bronski and Johnson 
[6], it seems remarkable that we are able to obtain the Whitham system ()4.12j) explicitly 
in terms of moments of the underlying periodic wave and the traveling wave parameters. 
Because of this, the verification that the hyperbolicity of the Whitham system correctly de- 
scribes the spectral stability near the origin proved to be much simpler and straight forward 
than cases previously considered by Oh and Zumbrun [20] and Serre. 

This issue of a missing equation comes up when there are more periodic solutions than 
would generically expected [26] , as arises when there is an integral of motion for the traveling 
wave profile ODE. Serre states on p. 262 of [26] that "a supplementary integral has a 
counterpart at the profile level", that is, such ODE integrals of motion often come from 
conserved quantities The "supplementary integral" he refers to is exactly an additional 
conservation law, which he proposes to use as above to fill the same missing equation it 
causes. That is, this same basic approach, by the argument of Serre, should work in a wide 
range of situations. 

In situations such as the KdV, where there may be still further conservation laws, 
there is flexibility in deriving the Whitham system, but all such derived systems must 
be equivalent. This would appear to correspond to existence for the Whitham system 
of additional conservation laws, or hyperbolic "entropies", another interesting possibility 
pointed out by Serre. 

We here just briefly investigate further the statement of Serre that supplementary PDE 



REFERENCES 



17 



integrals should imply also supplementary ODE integrals. We find justification for this 
statement in two very general situations. 

1. (more general). If by a supplementary integral, we mean an "entropy" in the con- 
servation law sense, i.e., a function rj(u) of the unknown and no derivatives, such that 
r](u) = d x (...), then this is indeed always true. (Note: the quantity J r]{u(x))dx is then 
conserved in time.) 

Proof. Proof: we then have additional traveling wave relation —srj(u)+q(u)+C = L(u, u x , ...), 
where C is a constand of integration and L is an expression of the same order as the 
traveling- wave ODE. But, combined with existing equations, we can then eliminate the 
highest derivative to get an algebraic relation between some lower derivatives, i.e., and 
integral of motion of the traveling-wave profile ODE. □ 

2. In the Hamiltonian case, both properties can be inherited from an overarching La- 
grangian structure, as in the multi-symplectic form system studied by many authors; see for 
example [21 E]. This may tie in an interesting way to the variational approach of Witham. 

We cite as open problems the rigorous verification of the Witham approximation at 
the level of behavior, either in the small e limit as in [llj . or in large-time behavior as in 
[231 1231 l25j I21j . Moreover, it is expected that a similar homogenization procedure can be 
used to justify the Whitham equations for the generalized Benjamin-Bona-Mahony equation, 
where the low frequency asymptotics were recently computed in a similar way as to the 
gKdV case: see [T3] . 

Acknowledgement. Many thanks to Jared Bronski for providing several useful insights 
into the connection of this work to that presented in [28J by Whitham. Also, we thank the 
reviewer for a careful and detailed check of the formulas and derivations presented in section 
four. 
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